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1.  Background 


This  note  presents  a cooputatlon  technique  for  efficiently  imple- 
menting a common  inventory  model  for  the  case  when  lead  time  demand  is 
assvuned  to  be  negative  binomial.  There  is  empirical  evidence  to  justify 
use  of  the  negative  binomial  for  m.  .ly  items. 

The  model,  which  is  described  In  the  next  section,  is  currently 
being  used  by  the  Department  of  Defense  in  wholesale  management  of 
secondary  Items  [4] . In  implementing  this  model  It  Is  necessary  to 
determine  expected  time-weighted  backorders;  to  determine  the  probability 
of  a stock-out;  and,  to  determine  the  optimum  reorder  warning  point  (R) 
for  a given  reorder  quantity  (Q). 

The  computation  technique  presented  uses  the  Camp-Paulson  approxima- 
tion to  the  negative  binomial  [i]  to  accomplish  these  objectives.  In  order 
to  apply  the  Camp-Paulson  approximations,  the  general  forms  of  the  in-r 
ventory  model  equations  must  be  transformed  using  identities  peculiar 
to  the  negative  binomial.  These  identities  are  proved  in  Appendix  1. 

2.  Model 

We  will  examine  an  Inventory  system  which  has  a central  stoekage 
point  from  whom  customers  requisition  stock.  The  inventory  system  is 
examined  using  the  following  assumptions: 

a.  As  demands  are  received  they  are  filled  from  available  stock. 

h.  When  stock  is  insufficient  to  fill  a requisition,  the  unfilled 
portion  is  backordered. 

c.  When  the  assets  (stock  on-hand  plus  on-order)  reach  the  reorder 
warning  point,  R,  a procurement  action  is  initiated  for  the  amount  Q.  The 
check  of  the  reorder  warning  point  is  done  on  a continuous  review  basis. 


^^^See  for  example  "Evaluation  of  Several  VSL/EOQ  Models,"  R.L.  DeeMr 
and  U.K.  Kruse,  AMCIRO,  ALKC,  Ft.  Leo,  Va.,  May  1974,  AD  78194d,  and 
reference  (3]  of  this  report. 
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d.  Procurements  will  raise  the  asset  level  to  the  requisition 
objective  (R+Q) . 

e.  Negative  reorder  points  are  not  allowed » i.e.,  R^o,  and  the 
Item  will  be  procuredt  l.e.»  Q ^ 1. 

£.  The  procurement  lead  time  Is  deterministic. 

g.  Demands  during  the  procurement  lead  time  follow  a negative 
binomial  probability  function. 

We  wish  to  minimize  the  annual  cost  of  holding  and  ordering  stock 
at  the  stockage  point  subject  to  a constraint  on  the  time^welghted 
requisitions  short.  The  holding  cost  Is  expressed  as 

H (R  + |) 


where  H Is  the  holding  cost*  e^^ressed  in  dollars  per  unit  per  year. 
There  will  be  D/Q  orders  per  year  and  hence  the  ordering  cost  la 


where  F Is  the  variable  coat  associated  with  Initiating  a procuramant 
action. 

Hadley  and  Whitln  [6]  give  the  general  axpresalon  for  expected  units 
backordered  as 


1 • 

B<R)  - i I (1-R)  lC(lflit)  - G(14Qfljt)l 
1-R 


(1) 


where  G(y;t)  It  tho  probability  of  having  y or  nora  demands  in  time  t. 
Expected  requisitions  short  era  approximated  as  the  expected  units  back- 
ordered divided  by  the  avarege  requisition  alssi  S.  [7] 

The  objective  of  the  model  la  to  ainialxe  total  annual  variable  coat* 
Kathematlcally,  this  will  be: 


minimise: 


IVC(R) 


TA  g 


(2) 
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vhere  X Is  an  implied  requisition  short  cost.  This  value»X,  can  be  deter- 
mined In  one  of  two  ways.  For  example » when  the  maximum  average  number  of 
requisitions  short  Is  specified  for  a group  of  items,  this  will  imply 
a shortage  parameter  and  this  affects  the  Inventory  funds;  or,  If  the 
total  funding  level  for  safety  level  investment  Is  fixed  for  a group  of 
Items,  this  iQq>lles  a shortage  parameter  and  establishes  the  performance 
that  can  be  achieved. 

The  optimal  reorder  quantity,  Q,  is  computed  in  some  non-optlmal 
manner.  For  example,  the  optimum  Q could  be  approximated  by  the  Wilson 
Q or  some  other  approximate  expression. 

Only  the  optimal  R will  be  solved  for  exactly.  Since  TVC  is  convex 
in  R,  difference  equations  must  be  solved  for  optimum  R.  Thus,  the 
optimality  condition  is 

TVC(R)  - IVC(R-l)  < o < TVC(R+1)  - TVC(R)  (3) 

3.  Fundamental  Identities 

The  general  form  of  the  negative  binomial  probability  function  is 

- r(r)  r<yfi)  <4) 

£<y;r)  ■ o y < o 

where  o^p^l,  q«l<-p,  r(xfl)  « x!  if  x is  an  integer  and  r(»fl) 

■ X r(x) . 


The  mean  and  variance  of  the  negative  binomial  distribution  are 
given  bys 


Evaluation  of  (3)  depends  on  the  follovlng  identities.  Let  i be  an 
Integer  and  R ^ o for  all  the  results.  Then* 

1 f(l;r)  - ^ f (l-l;r+l)  (5) 

* 

Z 1 f(ljr)  - ^ P(R-1;  r+1)  (6) 

1-0  P 

E 1^  f(i;r)  - ^ F(R“2;  r+2)  + F(R-1;  x+lA  (7) 

1-0  P P ' 

R 

E F(l;r)  - (R+1)  F(R;r)  - ^ F(R-lj  t+1)  (8) 

1-0  ** 

£ 1 F(i;r)  - P(R;r)  - ~ f(R-2;  r+2)  (9) 

1-0  ^ 2 p p 

where  F(y;r)  is  the  cumulative  distribution  function.  The  derivation  of 
these  identities  la  given  in  Appendix  1. 

Furtheraore«  F(ytr)  and  G(yit}  ate  related  by 

C(yfl{r)  - I - F(y;r)  (10) 

and 

M M**! 

E G(l;r)  - 1 E U-F(lsr)J  (11) 

1-0  1-0 

Probability  of  Stock  Availability 

The  evaluation  of  expretalon  (3)  requlret  the  calculation  of  piobability 
of  stock  availability  at  a function  of  R*  Xn  particular 

B(R-l)  •>  B(R)  - 1 - a (12) 

share  a denotea  probability  of  availability  of  stock  (sea  Appendix  2). 


2 

The  general  expression  for  the  probability  of  stock  availability, a, 
is  given  In  Hadley  and  Whltln,  [6}.  It  is 


♦ 


Transforming  this  to  the  cumulative  distribution,  F(y;r)  « 1 - G(y+l;r), 
and  using  the  Identities  of  section  3,  the  expression  for  the  availability 
of  stock  Is 

- ( R+Qf!  R+Q^-l  R-1  R-1  \ 

o-1-ijl  + Z 1-Z  7(iit)  - 1 - I 1 + I F(lir)  J 
^ V.  i»o  l»o  !■© 

• 1 - i {(R+Q)  - (R4Q)  F(R+Q-lir)  + 

^ F(R-HJ-2 {t+D  - R + R P(R-l;r)  - ^ F(R-2;  r+1)} 

■ 1 - ~|r  F(R-ltr)  - ^ F(R-2;  r+1)  + Q - 

(R+Q)  FCR+Q-ljr)  ^ F(R4<J-2{  t+1)^ 

* I {(R+Q)  F(R+Q-ljt)  - ^ F(R4<l-2i  r+l)  - 

R P(R-l;r>  + ^ F(R-2;  rfl)^  (14) 

5.  Expected  Humber  of  Backorders 

The  expression  for  the  expected  nuaiber  of  backorders  is: 


^Is  expression  is  generally  only  true  for  requisitions  of  unit  sire. 

An  expression  fer  non-unit  requisition  sire  was  developed  by  E.A.  Silver, 
Chung-ttei  Ho,  R.t.  Oeemor,  "Cost-Minittieing  Inventory  Control  of  Items 
Having  a Special  Type  of  Erratic  Demand  Fattam,"  IHFOR,  Vol  9,  Ho.  3, 
Movetober  1971. 
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B(R)  = r I (1-R)  [G(i+l;r)  - G(i+Qfl;r)l 
^ i»R 

Now  the  subtraction  within  the  brackets  is  accomplished  and 

® “ Q {ltG(R+2;r)  - G(R+Qf2;r)]  + 

2(G(R+3;r)  - G(R4<H3;r)l  + 

+ Q [G(R+Q<-l;r)  - G(R+2Qfl;r)]  + 

(0+1)  (G(R4<}+2jr)  - G(R+2Q+2;r)J 

+ ~ } 

■ I ■{c(R+2ir)  + 2G(IM-3it)  + — + 

QG(R+<}+l;r)  + QC(R-KH2;r)  + 

QG(R4Qf3;r)  + — J- 
-R+Q 


(15) 


- 4{s  (i-r)  G(i+ljr)  + Q i:  G(l;r)\ 

1"R+1  1-R-K)f2  / 

B • r (i-R)  U-F(itr)J  + I G(i;r)  - 

R+Qfl  1 

I C(l;t)J  f 
i-l  J 


1/R+O 

(i-R)  U-F(l;r)]  - Z (i-R)  [l-F(l{t)J 


•1-0 


R 

Z 

l»o 

ra  R'^O  ^ 

+ Ol“^  + l-  l-  E U*  P(l;r)J  V 
l«o  J 


The  mean  part  (rq/p)  of  this  oMpresslon  is  based  on  the  expected 
value  expression  found  in  (8J*  l.e., 


X « £ Pr(X  > 1}  and 
1-1 


X + 1 - E Fr(X  > 1). 

1-0 


Continuing  the  algebraic  operations  the  expressions  for  the 

expected  backorders  is 
* 

1 f R+Q  m 

B - X i[l-F(iir)l  - R Z [1  - ?(l;r)l  - 
^ *"i-o  i-o 

R R 

Z i{l-F(i}r)l  + R Z [1  - F(i;r)l  + 
i«o  i“0 

R4Q  N 

Q[^-S  tl-F(i;r)j} 

P loo 

J ^ F(R+<}-2;r+2)  - R(R4Q+1)  + 

R(R+Q+1)  F<a+Qir)  - F(R-HJ-l;r+l>  - 
Ri|til  ^ MBiil  p(R.,)  . |1  Xl^la  P(R-2{r+2)  + 

R(Rfl)  - R(RH)  F(R}r)  + R ^ P(R-l;r+l)  ♦ 

Q I ^ - (R4Q<-1)  + (R4<^fl)  P(R+<l;r)  - ^ P(R4<j-lir+l)  J 
B - ^ "^KR-M))^  + R+Ql  P(R4Q}r)  - RCR-l-l)  P(R;if)  + 

la  Xsiila  tFC\4q-2;r+2)  - P(R-2;t+2)J  + 

2 tRP(R-l;r+l)  - (R4<J)  r(R-fQ-l;x4).)]  + 

Q (2  (2R4<lfl)j  ^ (16) 
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6.  Reorder  Warning  Point  CompuT  atlon 


To  olninlze  the  total  annual  variable  cost  It  will  be  necessary  to 
have  a value  for  the  shortage  cost.  We  will  assume  that  one  of  the  two 
managerial  methods  described  In  Section  2 will  be  used  to  establish  a 
value  for  X. 

In  Section  2 we  also  asjtuaed  that  a value  for  Q would  be  given.  Hence, 
the  reorder  warning  point  is  the  only  variable  left  to  be  evaluated.  Let 
ua  proceed  to  solve  (3)  for  R. 

TVC(R)  - TVC(R-l)  i 0 < TVC(R+1)  - TVC(R) 

H(R  + f)  + - U(R-1+|)  - 

P|  -X  ^ o i U(R+1+|)  + P I + 

,mn.  .H(R+f).i^  -x^ 

«4kich  ylelda 

B(R-l)  - B(R)  i S—  < B(R>  - B(R+1)  (17) 

s 

the  optimal  R which  satisfies  (17)  will  be  the  largest  R such  that^ 

B(R-l)  - B(R)  < -4~  (18) 

- 8 . . ■ 

From  (12)  we  get 

1 - e l-f-  (19) 

S 

% ' ' 

See  (6]  for  derivation  of  this  statement. 
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(14)  gives  the  availability *0,  as  a function  of  R.  Since  R can 
not  be  solved  explicitly  we  have  to  search  to  find  the  largest  R such 
that  (19)  holds. 

Although  we  are  seeking  the  optimum  reorder  warning  point » RWP,  It 
will  be  easier  to  find  the  optimum  safety  level  wTiere  RWP  Is  made  up  of 
the  expected  demand  during  the  procurement  lead  time,  y.  , plus  a safety 
level,  SL,  to  provide  protection  for  the  random  nature  of  demand.  The 
RWP  can  be  written  as 

R • + SI. 

where  Uj^  can  be  considered  a constant  for  any  given  particular  calculation 
of  R and  the  SL  is  the  portion  of  the  RWP  which  will  vary  with  the  degree 
of  protection  desired. 

It  Is  possible  that  the  optimum  safety  level  will  turn  out  to  have 
a negative  value.  However,  most  managers  are  reluctant  to  implement 
negative  safety  levels.  Therefore,  for  practleal  considerations  a lower 
llalc  of  zero  will  be  established. 

Let  SL  be  the  optioum  safety  level  and  SL  be  any  non^nagatlva  value. 
Then, 

TC(SL)  I TC(SL*)  (20) 

where  TC  represents,  as  before,  the  total  annual  expected  coat  which  In- 
cludes holding,  ordering  and  penalty  costa.  Thus, 

* 0 n B(V|+SL  ) 

TC(SL)  >.  + SL  *1*  |)  + P|  + 1 Kg 

If  a positive  term  is  dropped  from  the  right  hand  side  of  tha  Inequality, 
tha  telatlonshop  still  holds,  i.e.. 


TC(SL)  >.  + SL*  + |>  + 


and 


B(p  + SL)  * 0 n 

+ SL  + J)  + + X ^ > H(Wj^  + SL  + f)  + 


which  yields 

B(v,+SL)  ^ 

H(Uj^  + SL)  + X ^ ^ + SL  ) 


Since  SL  is  evaluated  at  any  non-negative  value*  let  SL  •»  o.  Then, 


B(y.) 

^ S 


> + SL  ) 


and 


i-  > SL* 

H S - 


(21) 


This  expression  provides  an  upper  bound  for  the  safety  level. 

7.  Camp-Paulson  Approximation 

The  evaluation  of  the  upper  bound  of  the  safety  level  (21)  and 
the  inequality  which  must  be  satisfied  for  the  optimal  R,  (19),  require 
the  evaluation  of  the  negative  binomial  cumulative  probability.  For 
larger  values  of  R+Q  this  could  require  a significant  amount  of  time 
fcr  computation  purposes.  Two  probability  distributions  are  also  required; 
one  for  parameter  r,  and  one  for  parameter  r+1. 

Bartko  [1]  has  shown  that  a good  approximation  for  the  cumulative 
negative  binomial  is  a transformation  of  the  negative  binomial  variate 

4 

to  a standard  normal  variate  , i.e., 

NB(x;r,p)  - t (^J=)  (22) 


where 


4 

The  article  in  Technometrics  had  a typographical  error  which  has  been 
corrected  here. 


U 


y 


. 9 X + 8 
x+1 


[- 


1/3 


^(9r-l)  ^ p(yfl) 


-13 ■ 

p(yM) 


2/3 


x+1 


To  calculate  the  necessary  probabilities  for  the  optimum  R now  will 
simply  require  6 evaluations  to  get  the  maximum  safety  level  (equations 
21  and  16)  and  4 evaluations  for  each  R value  which  Is  appraised  In  the 
search  routine  (equations  19  and  14) . 
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APPENDIX  I 


IDENTITIES 

Section  3 gave  identities  which  were  instrumental  in  the  development 
of  the  solution  of  the  optimal  reorder  warning  point.  These  identities 
will  now  be  derived.^ 

The  negative  binomial  probability  function  is« 


£(y;r) 


-r(yty) 

r(r)  r(y+i) 


pV 


f(y;r)  - o 


y 

y < ‘o 


(23) 


The  mean  and  variance  of  the  negative  binomial  are* 


w - 


B. 

P 


3^  - B 


The  identities  are  developed  in  sequence  since  most  of  the  results 
depend  on  previously  derived  expressions.  First,  let  i be  an  integer. 
Then 


1 f(ijr)  - i 


r(r<»i) 
r(r)  rCi+1) 


r(r>U-i)i 


B 

p 


r(rfi) 

rr(r)r<i) 


p 


r+1  i-1 


P 


.jsmL 

r(r+i)r(i) 


„r4-l  1-1 
p q 


1 f(l;r)  - ^ £(i-l;rfl) 


(24) 


5 

Chung-Mci  Citem  did  much  of  the  early  work  in  the  development  of  these 
identities. 
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The  second  expression  is  Che  sum  of  the  first.  Thus, 

R R 

Z i f(i;r)  » L i f(l;r) 


R R R 

Z F(i;r)  - s I f(jjr) 

1-0  J-o  l«j 

R R 

- I f(j{r)  Z 1 

^ j-o  i»j 

R 

- Z f(j;r)  (R-j+1) 
j-o 

R R 

- (R+l)  Z f(j;r)  - Z j £(j;r) 

j-o  j-o 

R 

Z F(l;r)  - (R+l)  F(R;r)  - F(R-l;r+l)  (27) 

1-0  P 


The  last  Identity  of  concern  Is 

R R 1 

Z 1 F(l;r)  - Z 1 Z f(j;r) 

1-0  1-0  j-o 

By  changing  the  order  of  susuoatlon  ve  get 

R R R 

Z 1 F(l;r)  - Z 2 1 f(jjr) 

1-0  j-o  1-j 

R R 

- Z f(j;r)  Z 1 

J-o  1-j 

R rR  j-l 

■ S fCJsr)-}  Z 1 - Z 1 > 

J-o  t imi  J 

.1  ,(j„){asp  -.u^} 

£(j;r)iJ  (j*-j>  f(jir) 

^ j-o  * j-o 

• ^ - F(Rit)  - ~ F(R-2;r+2)  + 

r(R-lirfl)J  + I ^F(R-Xjt+X) 
t 1 F(l;r)  - F(R;t)  - I ^ ittyA...  r(R-2jw2)  (28) 

1-0  ^ 2 p p 

X5 
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APPENDIX  II 


PUTIONSHIP  OF  AVAILABTT.TTY  STOCK  AND  EXPECTED  BACKORDERS 

Consider  two  inventory  systems  which  differ  only  by  their  reorder 
warning  points,  one  using  R.  the  other  R-1.  Let  these  systems  start 
with  0 stock  at  time  t»o  and  subject  them  to  demands  starting  at  t«o 
to  t«»,  and  have  them  replenish  according  to  the  rules  stated  earlier 
(Section  2)  for  t>  o.  If  i^(t),  and  Ij^_^(t)  are  the  net  on  hand  values 
of  the  systems  with  reorder  points  R,  R»i  respectively,  then  i (t)  - 

+ 1 for  t > PLT  and  - l^^_^(t)  for  t < PLT  where  PLT  is 

the  procurement  lead  time.  Define  a stockout  period  as  the  time  from 
when  I^(t)  reaches  zero  from  a previous  positive  value  until  it  next 
returns  to  o.  Let  t^(R)  be  the  length  of  the  stockout  period.  Apart 
from  the  backorders  which  occur  for  t < PLT  (which  are  the  same  for  both 
systems)  no  backorders  occur  in  either  system  except  during  stockout 
periods.  Further  if  b^(R).  b^(R-l)  are  the  time  weighted  backorders 
during  t^(R)  for  systems  with  R.  R-l  respectively,  then 

• b^(R-l)  - t^(R)  * 1 unit 
and 

1 N(X) 

B(R)  ■ lira  b.(R) 

^ 1.1  1 


where  N(T)  • number  of  atockout  periods  daring  the  Interval 
(PLT.  PLT + T).  Then 

B(R-l)  - B(R)  ■ 11m  t.(R) 

TiMIII»S»1.1  * 

But  this  la  Just  the  proportion  of  tine  no  stock  is  available  or  l-o. 
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